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ABSTRACT 

The general relativistic precession rate of periastra in close-in exoplanets can be orders of magnitude 
larger than the magnitude of the same effect for Mercury. The realization that some of the close-in 
exoplanets have significant eccentricities raises the possibility that this precession might be detectable. 
We explore in this work the observability of the periastra precession using radial velocity and transit 
light curve observations. Our analysis is independent of the source of precession, which can also have 
significant contributions due to additional planets and tidal deformations. We find that precession 
of the periastra of the magnitude expected from general relativity can be detectable in timescales of 
< 10 years with current observational capabilities by measuring the change in the primary transit 
duration or in the time difference between primary and secondary transits. Radial velocity curves 
alone would be able to detect this precession for super-massive, close-in exoplanets orbiting inactive 
stars if they have ~ 100 datapoints at each of two epochs separated by ~ 20 years. We show that 
CO ' the contribution to the precession by tidal deformations may dominate the total precession in cases 

where the relativistic precession is detectable. Studies of transit durations with Kepler might need to 
take into account effects arising from the general relativistic and tidal induced precession of periastra 
for systems containing close-in, eccentric exoplanets. Such studies may be able to detect additional 
planets with masses comparable to that of Earth by detecting secular variations in the transit duration 
O ' induced by the changing longitude of periastron. 

Subject headings: celestial mechanics — planetary systems 
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1. INTRODUCTION 

Following the discovery of an extra-solar planet around the solar type star 51 Pegasi (Ma yor fc Quelozlll995l ) there has 
been rapid progress in the detection and characterization of extra-solar planetary systems. The very early discoveries 
CO ■ have shattered our view on planetary sys tems, as certain systems exhibited short periods (51 P eg), high eccentric ities 
(e.g. 70 Virginis b. lMarcv fc Butlerlll996h . and massive planetary companions (e.g. Tau Boo b. iButler et alJll997h . 
Interestingly, systems with all these properties combined (i.e. massive planets with short periods , small semi-major 
\Q ■ axes, high eccentricities) have been also discovered (e.g., HAT-P-2b, XO-3b; iBakos et alj 120071 Uohns-Krull et all 
' l2008[ ). The high eccentricities are somewhat surprising, as hot Jupiters with short periods are generally expected to 
be circularized in timescales shorter than the lifetime of the system if the param eter Q, inversely proportional to the 
■ plane t's tidal dissipation rate, is assumed to be similar to that inferred for Jupiter (|Goldreich fc Soterlll966l : lRasio et al.l 
1199(1 . 

By virtue of their small semi-major axes and high eccentricities, the longitude of periastron u> of some of the newly 
discovered systems are expected to precess due to General Relativistic (GR) effects at rates of degrees per century. 
This is orders of magnitude larger than the same effect observed in Mercury in our Solar System (43"/century), which 
offered one of the cornerstone tests of GR. Furthermore, the massive, close-in eccentric planets induce significant reflex 
motion of the host star, therefore enhancing the detectability of the precession directly via radial velocities. 

In this work we explore the observability of the precession of the longitude of periastron with the magnitude expected 
from GR in exoplanets using radial velocity and transit timing observations. We also consider in this work the periastra 
precession due to planetary pert urbers and tidal deformations, which can ha ve contributions comparable or greater 
than that of GR. Previous works (Miralda-Escudc 2002; He vl fc Gladmanll2007f ) have explored some aspects of the work 
presented here in th e context of us i ng tim ing observations to detect terrestrial mass planets. We refer the reader to 
independent work bv lPal fc Koc sis (2008) that also explores the measurable effects of the periastra precession induced 
by GR. 

2. EXPECTED PRECESSION OF PERIASTRA 

Before continuing let us fix our notation. In what follows a will denote the semi-major axis of the Keplerian orbit 
of the planet-star separation, e its eccentricity, P its period, uj its longitude of periastron, M* and i?* the mass and 
radius of the host star, respectively, and n = (GMtot/a 3 ) 1 ^ 2 is the Keplerian mean motion (orbital angular frequency), 
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TABLE 1 

GR PRECESSION OF EXOPLANETS WITH e > 0.1 AND GR PRECESSION RATES 

> 1° /CENTURY 



Name 


a 


e 




Period 


K 






TEP? 




(AU) 




(Mq) 


(days) 


(m sec -1 ) 


(deg) 


"/century 




HD49674 b 


0.058 


0.290 


1.060 


4.944 


13.7 


283.0 


1.576 





HD88133 b 


0.047 


0.133 


1.200 


3.416 


36.1 


349.0 


2.958 





GJ 436 b 


0.028 


0.159 


0.410 


2.644 


18.7 


339.0 


2.234 


1 


HD118203 b 


0.070 


0.309 


1.230 


6.133 


217.0 


155.7 


1.231 





HAT-P-2 b 


0.069 


0.507 


1.350 


5.633 


884.0 


184.6 


1.836 


1 


HD185269 b 


0.077 


0.296 


1.280 


6.838 


90.7 


172.0 


1.046 





XO-3 b 


0.048 


0.260 


1.410 


3.192 


1471.0 


-15.4 


3.886 


1 



where G is Newton's gravitational constant and M tot = M* + M p with M p the mass of the planet. The reflex motion 
of the host star is characterized in a similar manner. In this section we detail the expected mechanisms that will cause 
a precession in the value of w. 

2.1. General Relativistic Precession 

One of the most well-known consequences of General Relativity is that orbits in a Schwarzschild metric are no 
longer closed as is the case for the Kepler problem in Newtonian mechanics. The rate of precession of the longitude of 
periastron due to GR is given to leading post-Newtonian order by 

3GM* 

(see any general relativity textbook, e.g. lMisner et al.lll973l for a derivation). With a expressed in astronomical units, 
P in days and the mass of the host star M* in solar units, uqr is given in units of degrees per century by 

7.78 [ M+\ / a \-i / P V 1 ro . 

We use equation [2] to estimate the ex pected precession f or currently known exoplanets as listed in the online California 
& Carnegie Catalog of Exoplanets (jButler et al.l I2006L version Nov 7 2007). We list in Table [T] all exoplanets that 
have e > 0.1 and that have loqk > l°/century 5 . From left to right, the columns in Table Q] record the name of the 
exoplanet, its semi-major axis, its eccentricity, the mass M* of the host star, the period, the velocity amplitude 6 , the 
longitude of periastron u>, the calculated wqr and finally a flag that is 1 if the planet is transiting its host star and 
otherwise. The condition on e is in order to consider only systems where GR effects can be constrained with sufficient 
confidence. Additionally, we restrict ourselves to stars which have evidence currently for a single exoplanet in order 
to avoid the complicatio ns arising from the precession of u> induced by other planets (see below) 7 . W e have added to 
the exoplanets listed by iButler et alj (|2006f ) the recently discovered XO-3b (|Johns-Krull et aT]|2008f ). which has the 
largest predicted wqr of all exoplanets listed in Table [TJ 

As can be seen in Table [TJ seven systems have wqr > 1° /century, with 3 of them being transiting systems. In 
timescales of a few tens of years, the longitude of periastron of the systems is expected to shift in these systems 
by Sw ^ 0.5°, a change that, as we will show below, may produce detectable effects. GR effects are not the only 
mechanisms that can cause a shift in lu though, so we now turn our attention to additional mechanisms. 

2.2. Stellar Quadrupole, Tides and Additional Planets 

Bes ides the GR precession discussed above, u> can precess due to the presence of additional effects. Miralda-Escude 
(2002) discussed the observability, using the duration of transits and the period between transits, of changes in u> caused 
by a stellar quadrupole moment an d perturbations from ot her planets. These effects, some of which were discussed 
using more accurate calculations by lHevl fe Gladmanl (|2007h . may additionally cause a precession of the orbital plane. 
Additionally, tidal deformations induced on the star and the planet can also produce a secular change on u>, an effect 
not considered in the studies ment ioned above. The effect of apsidal motions induced by tidal deformations is a well 
known effect in eclipsing binaries (|Sternd[i939l ; iQuataert et al.1 1996: S meyers fc Willemal2001l) . and was included in 
the analysis of the planetary system around HP 83443 bv lWu fc Goldreichl (|2002h . Tidal deformations can produce a 
significant amount of precession in close-in exoplanets and should therefore be taken into account. 

5 The magnitude of all the effects discussed in what follows where uqh would manifest itself are incre asing functions of e. Systems 
with low e are therefore not releva nt from the point o f view of detecting GR effects. Hcvl & Gladman (2007) presents in their Figure 4 an 
estimate of ujqr f° r a H systems in Butler et al. (2006) and note the four systems with higher values. None of those are in Table [l] because 
they all have e < 0.03. 

6 This is often referred to as the semi-amplitude in the exoplanet literature. We choose to use simply amplitude in this paper in order 
to agree with the widespread usage in the physical sciences for the multiplicative factor in a simple harmonic oscillator. 

7 Adams & Laushlin (2006a b) studied the effects of secular interactions in multiple-planet systems including the effects of GR. They 
show that GR can have significant effects on secular perturbations for systems with favorable characteristics. 
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The precession caused by a stellar quadrupole moment is given to second order in e and first order in (i?*/a) 2 by 

3J2R 2 
2a 

where J2 is the quadrupole moment (|Murrav fe Dermottl [l999). In units of degree/century this expression reads 



Wquad ~ „ 2 Tl, (3) 



( —Y 1 (Ji\ 



^ - 0-17 ) (o^AuJ r/ceutury]. (4) 

It is clear from this equation that for values of J 2 < 10 -6 similar to that of the Sun (iPireaux fc Rozelot|[2003l ) the 
value of Wquad is smaller than the value of u> expected from GR (see also lMiralda-Escudell2002f ). We will therefore 
assume in what follows that w qua d is always negligible in comparison with ugr- 

The tidal deformations induced on the star and the planet by each other will lead to a change in the longitude 
of periastron which is given, under the approximation that the objects can instantaneously adjust their equilibrium 
shapes to the tidal force and considering up to second order harmonic distortions, by 

15/(e) (k2, s M p Rl k 2 , p M*R 5 p \ 
-tide « ^- [~^- + ) n, (5) 

where /(e) = (1 — e 2 )~ 5 [l + (3/2)e 2 + (l/8)e 4 ], and fc 2 . s , ^2, P are the apsidal motion constants fo r the star and planet 
respectively, which depend on the ma ss concentration of the tidally deformed bodies fStcrnc 193j|). For stars we expect 
^2,s ^ 0.01 (|Claret fc Gimene3 [l992'). while for giant pl anets we expect fc 2iP ~ 0.25 if we assume that their structure 
can be roughly described by a polytrope of index n « 1 (|Hubbardl [l984). Fo r the extreme case of a sp here of uniform 
mass density, the apsidal motion constant takes the value fc 2 = 0.75 (e.g., ISmevers fc Willemsll200il ). We see from 
Equation O that for close-in hot Jupiters the term containing fc 2!?) will dominate, and that the effect of tides on u> 
increases very rapidly with decreasing a. In units of degree/century equation [5] gives 

(JLY 1 (hA ( a r 5 (Rp\ B (MA (M± 



^j^j [W P ) [jfj < 6 > 

where we have introduced T = 1 + (R+/R p ) 5 (M p /M+) 2 (k2, s /k2, p ), which is 1 for the case of a close-in Jupiter. 
Assuming that k2, p ~ 0.1, e < 0.5, M p ~ Mj, M* ~ Mq, R+ ~ Rq and R p ~ Rj it follows that Wtide is of comparable 
magnitude as wgr- 

The precession of the periastra caused by a second planet, which we dub a "perturber", is given to first order in e 
and lowest order in (a/ 0,2) by 

3Af 2 a 3 
4M*a| ' 

([Murray fc DermottlH99l iMiralda-EscTidl 12002) , where M 2 is the mass of the second planet and a 2 the semi-major 
axis of its orbit. In terms of deg/century this expression reads 



^perturber ^^^^3^ V'J 



(P\ : M 3 (M£\ ( ii 2 



Wporturbcr -29.6^-j [-) {^) [7century]. (8) 

For a perturber with a 2 = 2a and and a mass similar to Earth orbiting a solar-mass star, we get that w pC rturbcr ~ 
3 x 10~ 7 n or w pC rturbcr ~ 0.7 deg/century for a P = 5 days planet. This is comparable to the precession expected 
from GR and therefore any detected precession of the longitude of periastron will be that of GR plus the possible 
addition of any perturber planet present in the system and the effects of tidal deformations (and generally a negligible 
contribution from the stellar quadrupole). 

In what follows we will discuss the observability of changes in lo using radial velocity and transit observations. As 
just shown, any precession is expected to arise by GR, the effect of additional planets or tidal deformations. The 
discussion that follows addresses the detectability of changes in uo independent of their origin. 

3. OBSERVABILITY OF PERIASTRA PRECESSION IN EXTRA-SOLAR PLANETS 

3.1. Radial Velocities 

The radial velocity of a star including the reflex motion due to a planetary component is given by 

v T (t) = v +K[cos(u + /(i-io)) + ecos(u;)], (9) 

where vq is the systemic velocity, to the time coordinate zeropoint, / the true anomaly and K is the velocity amplitude 
which is related to the orbital elements and the masses by 
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TABLE 2 

Results of Radial Velocity Curves Fit 
Simulations for = 100, o- bs = 2 

M/SEC AND (Tjittor = 4 m/sec 



K 


e 




"20 


(m sec -1 ) 




(deg) 


(dcg/century) 


100 


0.10 


3.77 


80.010 


100 


0.20 


1.85 


39.204 


100 


0.30 


1.29 


27.457 


100 


0.40 


1.05 


22.249 


100 


0.50 


0.95 


20.253 


100 


0.60 


0.86 


18.159 


1000 


0.10 


0.38 


8.059 


1000 


0.20 


0.19 


4.106 


1000 


0.30 


0.13 


2.768 


1000 


0.40 


0.10 


2.182 


1000 


0.50 


0.09 


1.959 


1000 


0.60 


0.08 


1.790 



where i is the orbit inclination. Fitting for the observed radial velocities of a star will give then direct estimates of 
vq, to, K, e, P and to. 

Our aim in this section is to determine if oj can be constrained tightly enough in timescales of tens of years or less 
in order to detect changes in u> of the magnitude produced by GR in those time-spans. In order to do this we have 
simulated data and then fit it with a model of the form given by equation [9] a total of 1000 times. We then recover 
the best-fit values of lu in all simulations and use that to estimate the probability distribution 4>(uj) expected under 
given assumptions. 

The systemic velocity vq and to are just zero-points that we set to in all our simulations, where we also set the 
time units such that P — 1 (note though that we do fit for all these quantities so that their effect on the fit propagates 
to the uncertainties of lu). By trying several values of to we have verified that the probability distributions recovered 
do not depend strongly on the particular value of uj, which we therefore fix for all simulations at an arbitrary value 
ujo = 135°. This leaves us with just two parameters to vary, namely e and K. 

Given e and K we generate N Q \, S data-points with times U uniformly distributed 8 within a period and then we add 
to each time a random number of periods between and 20. For transiting exoplanets the observations would have to 
be taken uniformly in time intervals excluding the transit. We then get the observed radial velocity from equation [9] 
as 



v r (ti) = K[cos(ujo + f(U)) + ecos(wo)] + G(0, a to t) (11) 

where 67(0, <7t t) is a random Gaussian deviate with mean and standard deviation cr to t- The latter quantity is obtained 
as erjf ot = (T^ hs + cr? itter , where cr bs is the random uncertainty for each measurement and critter is the noise arising from 
stellar jitter. Given the form of equation [9] a Fisher matrix analysis implies that the uncertainty in the longitude of 
periastron, satisfies the following scaling 

VuKN^aufiK- 1 . (12) 

This scaling allows us to perform a set of fiducial simulations for several values of e and use the results to scale to 
parameters relevant to a given situation of interest. We note that the simulations we performed verify that the scaling 
inferred from a Fisher matrix analysis is accurate. In order to measure the longitude of periastron oj with a reasonable 
degree of certainty the system clearly needs to have a significant amount of eccentricity. We restrict ourselves to 
systems with e > 0.1 and we simulate systems with e = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. 

For <7 bs we assume a typical high-precision measurement with o v,bs = 2 m sec" 1 . For the stellar jitter, we perform 
our fiducial simulations for a typical jitter of emitter = 4 m sec -1 ( Wr ight! [20051 : iButler et alj 120061 ). The limitations 



imposed by active stars and/or different precision on the radial velocity measurements on the recovery of wqr can be 
explored by using the scaling of with higher assumed values of emitter and/or <7 bs- Even though Equation[T2l renders 
multiple values of K redundant, we choose to present results for two values of K for illustrative purposes, namely 
K = 100 m sec" 1 and K = 1000 m sec" 1 . The former corre sponds roughly to Jupiter-mass exoplanets and is fairly 
representative of currently known systems (jButler et a l. 2006) 9 , while the latter corresponds to super-massive planets 
which as we will see are the class of systems which would allow the detection of ugr with radial velocities. Finally, we 
do our fiducial simulations for iV b s = 100, a value not atypical for well-sampled radial velocity curves available today. 

The results of the simulations are summarized in Table [2] From left to right, the columns in this Table record the 
assumed radial velocity amplitude K, the system eccentricity e, the expected uncertainty in the longitude of periastron 

8 We ignore in our simulations the Rossiter-McLaughlin effect for the case of transiting planets (Rossitcr 1924; McLaughlin 1924; 
Qucl oz et alj|2000h. 

9 See also http://www.exoplanet.eu 
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Fig. 1. — Distribution tj> of recovered angles of periastron u> — u>o, where ljq is the input angle, for the simulation with N obs 
K = 1000 m sec -1 and <r, = 4 m sec -1 . The distributions are shown for eccentricities e = 0.1,0.3,0.5. 



100, 



cr w for the simulated system, and finally a2o , which we define to be the value of lu that would be necessary in order to 
achieve a 3cr detection in the simulated systems in a time-span of 20 years when measuring u> in two epochs, each having 
Aobs observations. We note that we have performed Shapiro- Wilk normality tests in the recovered u> distributions. 
We found that all of them are consistent with normality and thus we are justified in using the dispersion a w to derive 
confidence levels. 

As an example, we show in FigureUJthe distributions of recovered u> for the case with N \, s = 100, K = 1000 m sec -1 
and a j = 4 m sec" 1 for eccentricities e = 0.1, 0.3, 0.5. In the upper panel of Figure[2]we show the radial velocity curve 
for a system with K — 1000 m sec" 1 , e = 0.5 and u> = 135°, while in the lower panel we show the difference between 
that curve and a radial velocity curve having to — 136°. The simulations presented in Table[2]show that this difference 
can be detected at the 10-cr level (a^ ~ 0.1°) when using AT b s = 100 observations, each having an uncertainty of ~ 2 
m sec -1 . The high level of significance can be achieved in this case thanks to the assumed low-level of jitter and the 
super-massive nature of a system with that semi- velocity amplitude. 

In order to validate our simulations against uncertainty estimates obtained with real data, we have run a simulation 
with Af bs = 20, e = 0.517, K = 1011 m sec -1 , u> = 179.3° and (Tj = 60 m sec -1 . These obse r vation al conditions and 
orbital parameters are appropriate for the observat i ons of HAT-P-2 b reported bv lBakos et al.l (|2007f ). Our simulations 
for this case return a w — 3.8°, while iBakos et al.l ()2007l ) report Co>HAT-P-2b = 179.3 ± 3.6°, in very good agreement 
with our estimate. We conclude that our simulations return realistic estimates of the uncertainties in ui. 

The simulations presented in Table [2] show that in some cases the precession of periastron detectable in 20 years, ct2o, 
is comparable to the values u>gr of currently known systems listed in Table [TJ For close-in, eccentric, super-massive 
planets (K ~ 1000 m sec -1 ) about 100 observations per epoch are sufficient, while for Jupiter-mass exoplanets 
(K ~ 100 m sec -1 ) an unrealistically large number of radial velocity observations, A Q b s = 10000 per epoch, would be 
needed. Therefore, observations of different epochs of radial velocities with a time-span of ~ 20 years would detect 
the variations in the precession of exoplanet periastra induced by GR in some currently known super-masive systems. 
The estimates above assume a typical level of stellar jitter. For active stars the number of observations have to be 
increased in proportion to the dispersion characterizing the jitter (see Equation I12p . Stellar activity is therefore an 
important limitation to detect changes in u) using radial velocities. All in all, radial velocity studies of exoplanets will 
be generally able to ignore the effects of GR prece ssion as they will usual l y be w ell below a detectable level. 

We note in closing that iMiralda-Escude! (|2002f ) and iHevl fc Gladmanl (|2007[ ) consider using radial velocities and 
transit timing observations in order to measure the small difference between the period observed in the radial velocities 
and the period between primary transits. Both works conclude that this is not a competitive method and so we will 
not consider it further here and refer the reader to those works for details. 

3.2. Duration of Transits 

As the planetary orbit acquires a significant eccentricity e, the duration of the primary and secondary transits are 
no longer equal and acquire a dependence on the longitude of periastron of the system . An explicit expression for 
the duration of a transit D in the eccentric case was derived bv lTinglev fc Sackettl (|2005l . their equation 7) under the 
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Orbital Phase 

Fig. 2. — (Top) Radial velocity curve as a function of orbital phase for a system with K = 1000 m sec -1 , e = 0.5 and u> = 135°. 
(Bottom) Difference between the curve in the top panel and a radial velocity curve that is identical except in that it has ui = 136°. Note 
the different y-axis scales used in the different panels. 

assumption that the distance between the planet and the star does not change significantly during transit. It is given 

by 

1/3 



D = 2Z{R+ + R p )- 



(l + ecos(/ t )) \2nGM tot 



P 



(13) 



where 



1- 



ri cos^ i 
(i?* + R p ) 



= v 7 ! -b 2 



(14) 



is a geometrical factor related to the impact parameter b. i?* and R p are the radii of the star and planet respectively, 
M tot = + M p and r t and f t are the radii and true anomaly at the time of transit, respectively. The latter clearly 

depends linearly on uj and we therefore have ft — &>■ The logarithmic derivative of the duration of a transit is 



dhiD/dt 



wesin(/ t ) 



1 



b 2 



l-b 2 



(15) 



(l + ecos(/ t )) 

In Figure[3]we show the quantity (1 /&))dln D / dt for a central transit (6 = 0) and for e = 0.1, 0.3, 0.5, 0.7, with higher 
e giving higher values of \(1 / io)d\n D / dt\. The change in the duration of an eclipse for a small change in u> given by 
uiSt is simply SD ~ DcuSt [6j~ 1 dh\D /dt]. For u>St ~ 0.5 x 10~ 2 rad, appropriate for the expected change in oj for 
wgr ~ 3 deg/century over 10 years, we have that 6D ~ 0.075Z? x 10~ 2 for e ~ 0.3, which translates into 8D ~ 10 sec 
for a transit duration of D ~ 0.15 day which is typical for the sytems that we explore in this work. 

A very interesting feature of Equation 1151 is its dependence on the impact parameter b. First, d\nD/dt vanishes 
for b = 1/V2- This behavior is possible due to two competing effects which cancel out exactly for that value of 
b: an increase/decrease in the star-planet separation causes both and decrease/increase in the path-length of the 
planet across the disk of the star and a corresponding decrease/increase on the velocity across it, which implies a 
slower/faster crossing-time. Secondly, for systems with values of b close to 1 (i.e., near-grazing systems), the value of 
dhiD/dt increases greatly. We will come back to near-grazing systems below; in what immediately follows we will 
quantify the accuracy to which we can determine the transit duration D. 

If the times of beginning of ingress and end of egress are denoted by ti and t e respectively, then the duration of 
a transit is given by D = t e — ti and the uncertainty in the duration is a 2 D — of. + of (assuming no correlation 
between ti and t e ). If the ingress has a duration of Ai, (i.e. the time from first contact to second contact) then a linear 
approximation of the flux during this time can be written as F(t) = Fq(1 — it — U)(Rp / R±) 2 / Aij) , where Fq is the 
out-of-transit stellar flux and we ignore the effects of limb darkening. If the light curve is being sampled at a rate T per 
unit time then we have N = AtiT photometric measurements during the egress. If each photometric measurement has 
a fractional precision cr p h, and assuming R p and i?* are known, a least-squares fit to the photometric series will allow 
the determination of ti with an uncertainty cr ti = cr p hAti(i?*/i? p ) 2 A^ -1 / 2 = tT p h(A^/r) 1 / 2 (i?*/i?p) 2 . The uncertainty 
in the time of egress is given by a similar expression replacing Ati by At e . 
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Fig. 3. — (l/ui)dln D /dt as a function of the true anomaly at the time of transit ft for a central transit (p = 0). The different curves 
are for difference eccentricities e = 0.1,0.3,0.5,0.7, with higher e giving higher values of \(l/u)d In D/dt\ . The extrema in this figure are 
at values of the true anomaly at the time of transit of ft = ±arccos(— e). 



If we assume that the ingre ss and egress tim es are equal 10 , 
get that for a central transit (jFord et al.ll2008l ): 



and insert explicit expressions for Atj w At e and D, we 



op 
D 



40.7cr ph yi + ecos/t 

(!_ e 2)l/4 



N tr T \R 



R„ 
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where N tr is the number of transits observed, T is to be expressed in units of min 1 , and we have neglected factors 
including (1 + fi), where /i = M p /M+. We note for later reference that the uncertainty in the central transit time 

t c = 0.5(< e + U) is cr c « a D j 

Th e highest photometric precision in a transit light curve has been achieved with HST ([1 rown et al.ll2001t|Pont et all 
l2007h with (Jph ~ 10 4 . For parameter s appropriate to HP 209458 Equation [TBI gives ao ~ 5 sec. The uncertainties 
in the central time er c ~ <td reported in lBrown et al. (2001) arc of the same order and we therefore deem Equation [TBI 
to be a reasonable estimate of the expected uncertainties 11 . We will take this to be the highest precision currently 
possible for facilities su ch as HST where on l y a few transits a re typically observed. 

The Kepler mission (jBorucki et alj [20 03; Basri et ahl 120051 ) will observe ~100,000 stars continuously with cr p h ~ 
4 x 10~ 4 for a one minute exposure of V = 12 solar-like star and expects to find a large number of close-in "hot 
Jupiters" . The precision attainable in one year for a 1-minute sampling of a Jupiter-mass system with P = 5 days 
orbiting a V — 12 solar-type star is (Jd/D ~ 1.1 x 10~ 4 or cd ~ 1.5 sec for a 0.15 days transit duration. 

Kepler should therefore be capable of detecting changes in the transit duration D due to GR within its mission for 
some eccentric, close-in systems, and certainly when coupled to follow-up determinations of D with a facility delivering 
a precision similar to what HST can achieve. Of particular interest will be systems that are near-grazing, due to the 
factor of (1 — 6 2 ) -1 in Equation [T5l If eccentric, close-in, near-grazing systems are found in the Kepler CCDs, they 
will be subject to significant changes in D which should be observable. For example, a system with e = 0.4, b = 0.85, 
ft = 7r/2 and u>gr ~ 3 deg/century would have a change in D of ~ 9 sec during Kepler's lifetime for D — 0.075 days, 
a change which would be detectable. Of course, as b approaches one, D will tend to zero and the assumptions leading 
to Equation [TBJ will break down and make the estimated uncertainty optimistic, both effects which will counter the 
corresponding increase of dhiD/dt. 

3.3. Period Between Transits 

As already noted by iMiralda-Es cudc (2002|) and lHevl fc Gladmanl (|2007f ). as the longitude of periastron changes, the 
period between transits Pt will change as well. Periods are the quantities that are measured with the greatest precision, 



10 The ingress and egress times are not equal in general for eccentric systems, see Equation 7 in I Ford et al.l 1120081 . This small difference 
has no significant effect on the uncertainty estimates dealt with here. 

11 We note that Brown et al. (2001) warn about the presence of systematic effects not accounted for in the Poisson uncertainty estimate 
that could be of the same order as for their observations. 
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usually with uncertainties on the order of seconds from ground-based observations (e.g., the average uncertainty for 
HATNet planets discovered to date is 4.5 seconds). 
To first order in e the derivative of the transit period is given by 

P t = 4ne(^-\ sin(Mt) (17) 

where Mt is the mean anomaly at transit (jMiralda-Esc udc 2003) and is related to the true anomaly at transit ft to 
first order in e by M t — ft ~ 2esin/ f . For e = 0.1, a 5 day period, and wqr = 3 deg/century, the root mean square 
value of dPt/dt over all possible M t values is ~ 10~ 12 , which translates into a period change of ~ 2 x 1CP 4 sec in 10 

years. The period can be determined to a precision ~ a c N tr 3 ^ 2 or aoN tr 1 2^ 1 / 2 using the expression for ctd above 12 . 
Assuming the parameters for Kepler as above (V = 12 star, P = 5 day period) the precision achievable during 1 year 
is ~ 0.013 sec. 

Based on the numbers above we conclude that measuring significant changes in the transiting peri od in < 10 years 
timescales is not feasible. Our con clusions are in broad agreement with the analysis presented in iMiralda- Escudel 
(2002| and lHevl fc Gladmanl (|2007f ). who conclude that thousands of transits need to be observed with high precision 
in order to detect significant variations in P t . As there is no existing or planned facility that will allow to observe 
this amount of transits with the required photometric precision we conclude that measurements of P t will not be 
significantly affected by changes in id of the magnitudes expected to arise from GR or from the secular changes due to 
a perturber. 

3.4. Time Between Primary and Secondary Transit 

In the case where the exoplanet is transiting it may be possible to observe not only the primary transit, i.e. the transit 
where the exoplanet obscures the host star, but also the occulta tion, when the host star blocks thermal emission and 
reflected light from the exoplanet (e.g.. lCharbonneau et al.ll2005l ). If the time of the primary eclipse is given by t\ and 
that of the secondary by t2, the time difference between the two as compared to half a period P, At = t<i — t\ — 0.5P, 
depends mostly on the eccentricity e and the angle of periastron uj. Indeed, an accurate expression that neglects terms 
proportional to cot 2 i where i is the inclination angle and is therefore exact for central transits, is given by 




P I ecos(wWl — e 2 , ~ , , 

At = - n \ — + arctan n g 18 

(|Sterndll940[ ). Combined with radial velocities this time difference offers an additional constrain on e and u>, and in 
principle a measurement of At combined with a measurement of the differe nce in the duration of the secondary and 
primary eclipses can be used to solve for e and u) directly (see discussion in lCharbonneau et al.ll200 5). We note that 
iHevl fc Gladm an (2007, their §4.2) also consider secondary transit timings as a means to measure changes in u>. While 
their discussion is based on first order expansions in e instead of using the exact expression above and is phrased in 
different terms, it is based ultimately on the same measurable quantity we discuss here. 13 

Equation [TH] does not include light travel time contributions, i.e. it neglects the time it takes for light to travel 
accross the system. This time is given for a central transit by 

2a(l-e 2 ) , , 

c[l - {e cos ft) 2 } 

where f t is the true anomaly at the time of primary transit. In this section we will be interested in changes in At 
due to changes in ui. It is easy to see from the expressions above that dA ti LT I dto <C dAt/doj (ignoring the points 
where they are both zero). Therefore, changes in the time difference between primary and secondary transits will be 
dominated by changes in At and we can safely ignore light travel time effects in what follows. 

The precession of periastra caused by GR will change u systematically, while leaving e unchanged. Therefore At 
has the potential of offering a direct measurement of changes in oj. We use Equation [18] to calculate (l/P)dAt/du> as 
a function of lu. The result is shown in in Figure 2] for e = 0.1,0.3,0.5,0.7, with higher e yielding larger extrema of 
\{\/P)dAt/duj\. 

If lj ss 7r/2 or « 3-7I-/2, the changes At in a system with P = 4 days for 0.2 degrees of precession would be large, 
on the or der of 8 min for e ~ 0.5. Published observations of secondary eclip ses with Spitzer constrain t% to within 
~ 80 sec (jDeming et al.ll2006t iHarrington et alJl2007t ICharbonneau et al.ll2008l ). Uncertainties in ti will dominate the 

uncertainty in A t , so we assume <7a, ~ &t 2 ~ 80iV tr 3 ^ 2 sec, where N tr is the number of secondary transits observed. 

Determing At at two epochs separated by 10 years and observing N tr ~ 1 at each epoch with a facility like Spitzer 
would detect the changes in At due to GR with a significance > 4er in an eccentric system that precesses > 0.2 

12 The scaling N~ 3 ' 2 follows from describing t he central tr ansit times as tj = to + Pi, with i = 1, . . . , N tr and determining P using a 
X 2 fi t. The variance in the d erived P is (see, e.g.. IGouldll200l) a 2 p « (3/iV 3 )o- 2 + OiN' 4 -). 

13 Hcyl & Gladman ( 2007) consider the time difference between successive primary transits (Ato) and successive secondary (At^ ) trans its, 
following their notation. This difference can be expressed as Ato — At w ~ —{d&t/duj)5ui, where At is the quantity defined in Equation 1181 
and 8lo is the change in uj in one orbit. 
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Fig. 4. — (1/ P)dAt/dhj as a function of lo. The different curves are for different eccentricities e = 0.1,0.3,0.5,0.7, with higher e giving 
larger extrema of (1/ P)dAt/duj\. 

deg in a decade if they have u> around tt/2 or 3tt/2. While Spitzer is now entering the end of its cryogenic lifetime, 
future facilities might be able to ac hieve similar precisio n. Moreover, as new suitable systems are discovered with very 
hot atm ospheres (pM class planets. | Fortnev et al.ll2007l ) observations of secondary transits might be feasible from the 
ground (|Lopez-Morales fe Seagerll2007f ). 

Finally, we note that since ft = ±7r/2 when lu = or tt, measuring the change in At complements measurements of 
changes in primary transit duration D, in the sense that when one effect is not operating the other one generally is. 

4. ASSESSING THE DIFFERENT CONTRIBUTIONS TO lu 

If a change in lu is detected using any of the methods described above it would be interesting to determine its likely 
cause. As mentioned in § [2j the contribution to lu from a second planet or tidal deformations can be of comparable 
magnitude to that caused by GR, making the origin of a potential detection of a change in lu unclear unless the 
presence of a perturber and the magnitude of tidal effects can be assessed. In this section we address three points. 
First, we discuss how may one assess if part of a detected change in lu comes from a perturber. Secondly, we show that 
the a terrestrial mass (M® < M < 10) perturber may be detectable through the induced secular changes in lu over the 
value expected from GR and tidal deformations by using transit and radial velocity observations. Finally, we assess 
the relative contributions of GR and tidal deformations for the case where perturbers do not contribute significantly 

to LU. 

4.1. Constraining Potential Contributions to lu from w per turber- 

First, we note th at using a precise ra dial velocity curve with N a h s ~ 100 we can probe the presence of companions 
with M > 15M e (jNaravan et all [2005). As shown in § lower mass companions can still contribute significantly 
to lu, so it would be useful to have further information about possible companions in order to better interpret a 
measured lu. The presence of a second exoplanet in the system will not only cause secular changes in lu as discussed 
in jj (equation [71 ), but it may cause significant trans it-time variations even for perturbers with masses comparable to 
the earth (jAgol et al.ll2005t iHolman fe M urray 2005T). The transit-time variations will be on the order of seconds to 
minutes, depending on the orbital parameters and mass of the perturber. Constant monitoring of transiting exoplanets 
such as will be performed by Kepler will allow to identify systems that show significant transit time variations. 

We note that the presence of a perturber will n ot only change lu but it may also affect other orbital elements 
such as the orbit inclination (jMiralda-Esc udc 2002). These small changes in inclination have indeed been proposed 
as a met hod of detecting t errestrial mass companions in near-grazing tr ansit systems which are esp ecially sensitive 
to them (|Ribas et alj [20081. Following the formalism presented in §6 of iMurrav fc DermottJ (|1999f ) we expect that 
secular changes to e by a perturber satisfy \de/dt\ < (3/16)n(a/a 2 ) 3 (M 2 /M+)e 2 , where a 2 ,e2 and M 2 are the orbital 
semi-major axies, eccentricity and mass of the perturber respectively. For e 2 = 0.1, a 2 = 2a and Af 2 = 3 x 10~ 6 M* 
we have that de/dt <4x 10 -6 yr _1 for a system with P = 4 days. This is too small for changes in e to be detected 
with radial velocities in scales of tens of years and so we conclude that secular changes in eccentricities induced by 
perturbers will in general not be useful to infer the presence of perturbers. 
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Summarizing, the best way to try to constrain contributions of cj per turber to a detected change in to is to consider 
transiting systems posessing extensive photometric monitoring, allowing to probe the existence of transit time varia- 
tions. For near-grazing systems, the same photometric monitoring would additionally allow to probe for small changes 
in % due to a perturber. 

4.2. Using Secular Variations in lj to Detect Terrestrial Mass Planets. 

A very interesting possibility raised by the secular variation in oj induced by a "pert urber" is to use these va riations in 
order to infer t he presence of terr e strial mass planets. This possibility was studied by Miralda-Escudc (2002) and then 
followed-up bv lHevl fc Gladmanl (|2007f ). Secular variations in u, especially through their effect on transit durations 
(§ 13.21) . can offer an interesting complement to transit time variations as a means of detecting terrestrial mass planets 
in the upcoming Kepler mission. 

As we have seen in § [2j the precession due to GR and tidal deformations can be of comparable magnitude to that 
induced by a terrestrial mass perturber. It is germane to ask then to what extent will the uncertainty in the expected 
value of wgr and Wtidc limit the detectability of a perturber. 

We start by considering wqr ■ The fractional uncertainty in wgr is given by (er(tjGR )/<^gr ) 2 = {o-{M+) / M+) 2 + 
(a(a)/a) 2 + (a{P)/P) 2 + 4e 2 (<r(e)/(l — e 2 )) 2 , where we have ignored correlations. The fractional uncertainty in P 
is generally negligible, while the other quantities can be typically of the order of a few percent. We will therefore 
assume that ct(wgr )/wqr ~ 10%. It follows that to detect an excess precession caused by a perturber we need at 
least that w pe rturber — d>GR ^ 0.3c1)gr • It is easy to see that if cL> per turbcr ^ wgr and cj pcrtU rbcr is itself detectable, i.e. 
w pertur ber > 3(7^, where is the uncertainty in the measured w, then the uncertainty in the expected wgr will not 
spoil the significance of the detection. Using the equations presented in §[2] we get that 

Wperturbor _ 3.8 / Mq \ 2 ( d\ ( a \ ( M 2 



w G R (l-e 2 )\AhJ \a 2 J V0.05AUV V M ffi7 ' 

It is clear from this expression that for typical values of M* ~ M , a ~ 0.05 AU, e < 0.5 there will be values of a 2 
for which we have w pex turber > ^>gr an d for which the uncertainty in the expected precession from GR will not be a 
limiting factor in detecting terrestrial mass perturbers. 

We consider now wtide- The fractional uncertainty considering all parameters excepting k 2 , P and T and ignoring 
correlations is (a^tidcV^tidc) 2 = (a{P)/P) 2 + (5a(R p /a)/(R p /a)) 2 + (a{M p /M*)/ \M P /M*)) 2 + (/'(e)a(e)//(e)) 2 . 
We have expressed this in terms of the ratios M p /M* and R p /a as these quantities are more robustly determined 
observationally. As was the case above, the uncertainties in the variables can be typically a few percent and we 
therefore assume that cr(<j t ide)/wtidc ~ 10%. Using the equations presented in §[5] we get that 

^perturber _ 1-85 / Mq \ 2 f a \ 3 / \ 5 / Rj\ 5 / M p \ / M 2 \ 

Lome Tf(e)k 2 , p \AlJ \a 2 ) Vo.05AuJ \R P ) \ Mj ) \ M e / ' 1 ' 

In the case of Wtide we cannot measure the apsidal motion constant k 2 . As discussed in § [21 the value of k 2 for a giant 
planet is expected to be close to the extreme value of a uniform sphere, so we can conservatively assume k 2p = 0.75, 
a value that maximizes the expected Wtide- For a close-in Jupiter we have T ~ 1. Using these values and following 
the same reasoning as for wqr , it is clear from the expression above that for typical values of M* ~ Mq, a ~ 0.05 
AU, e < 0.5, R p ~ Rj, M p ~ Mj, there will be values of a 2 for which we have Wperturber > ^tidc and for which 
the uncertainty in the expected precession arising from tidal deformations will not be a limiting factor in detecting 
terrestrial mass perturbers. Note that due to the strong dependence on a, tides may become a limiting factor for very 
close-in systems. 

We consider now the detectability of w p0 rturbor with Kepler using the change in the transit duration D. As discussed 
in § 13.21 a 1-minute sampling of a Jupiter-mass system with Kepler for a P = 5 days planet orbiting a V = 12 
solar-type star will achieve a precision of an ~ 1-5 sec for a 0.15 days transit duration. We therefore set a change 
of 6 sec over 4 years to constitute a detectable 5D for this system, which translates into a detectable Co of ^detect = 
6/(DAt(il)~ 1 d\nD/dt)), where At — 4 years and we assume a central transit (see Equation [T5"|). 

We show in Figure [5] the expected value of cL> per turber as a function of a 2 /a assuming a = 0.05 AU, M* = Mq, and 
ft = 0.5ir as solid lines, one for M 2 = Mq and another for M 2 = IOMq. The dashed line marks the values of W(j etect 
for e = 0.4. The dotted lines mark the values of 0.3wgr for the same eccentricity, while the dash-dot-dot line marks 
0.3wtidc This figure shows that over a range of values of (a 2 /a) Kepler may be able to detect the presence of additional 
super-Earths (Mq < M < IOMq) if the primary planet has a favorable impact parameter and true anomaly at the 
time of transit. 

Finally, we consider the detectability of w per turber using radial velocities alone. We assume the primary planet is a 
super-massive system of M p — 10Mj, with the rest of the star and orbital parameters as assumed above. The results of 
§ l3.1l show that with 100 radial velocity observations at each epoch a precession of ~ 8 deg/century could be detected 
(at 3cr) in 5 years if the system has e = 0.4 and orbits an inactive star. We mark the detectable level in this case with 
radial velocities as dot-dashed line in Figure [5) From this Figure we can see that if enough observations arc in hand, 
perturbing planets with masses ~ IOMq may be detectable with radial velocities alone for the case of super-massive 
systems orbiting inactive stars. 

We conclude that measuring the precession of u) caused by a perturber above the value predicted by GR and tidal 
effects may lead to detection of additional planets with both transit and radial velocity observations (for transits, 
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(a 2 /a) 

Fig. 5. — The solid lines show the expected value of tj pe rturber as a function of 0,2/0, for M2 = Mq and M2 = lOMfg, assuming a = 0.05 
AU, M* = Mq and /t = 0.57T. The dashed line marks the lower value of io that is detectable with Kepler for e = 0.4 via the changes in 
transit duration. The dotted line marks the value of 0.3(Ijgr f° r the same eccentricity, while the dash-dot-dot line marks 0.3<I> t id o . Finally, 
the dot-dashed line marks the lower values of lu that would be detectable with 100 radial velocities at each of two epochs separated by 5 
years, for a perturber to a super-massive Mi = 10Mj primary planet orbiting an inactive solar-mass star (see text for more details). 

see also iMiralda-Escudl (|2002D ; iHevl fc Gladmanl (|2007t )V The examples given above are illustrative only; a detailed 
analysis is beyond the scope of this work. Of special interest will be to estimate the yield of terrestrial-mass planets 
expected from Kepler by measuring secular changes in transit durations. 

4.3. The Relative Magnitude o/wgr and ujtidc- 

If there are no perturbers causing significant changes in u>, these changes will be caused by a combination of GR and 
tidal deformations. We now briefly consider the relative magnitudes of these effects. Using the equations presented in 
§ [2] we find that 

cOgr_ = 4.8 / OT \ / a \4 / RA 5 ( Mp 

d/tkie T(l-e 2 )f{e) \k 2 , P J \0.Q5AV J \R P J \Mj 

For expected values of k2 :P ~ 0.25 and assuming R p ~ Rj, M p ~ Mj, a ~ 0.05 and e < 0.5 we see that the values of 
cL>gr and Wtidc are comparable. Note that there is a very strong dependence on a. For smaller values of a, a regime 
where the precession due to GR would offer better chances of being observable for eccentric systems, the contribution 
from tidal deformations may quickly become the dominant source of precession. The fact that the tidal contributions 
to the precession are expected to be generally significant limits the ability to directly extract the precession due to 
GR. The magnitude of the precession caused by tidal deformations is uncertain due to the need to know the planetary 
structure. Given that the expected precession due to GR is unambigous, in systems where the precession due to tidal 
deformations is expected to be dominant and detectable, we might be able to learn about the tidally driven precession 
of uj by subtracting the effects of GR from a measured d>. 

5. ANALYSIS OF SPECIFIC SUPER-MASSIVE SYSTEMS 

As illustration of the material discussed above, we analyze now in some detail two currently known systems that 
due to their large measured velocity amplitudes have the potential to have a me asurable to using all the techniques 
presented above: HAT-P-2 b (|Bakos et al.ll2007t lLoeillet et alj |2008) and XO-3 b (|Johns-Krull et al.ll2008t ). 

5.1. HAT-P-2 b 

HAT-P-2 b is an especially interesting system due to its very high eccentricity, high mass of « 9Mj and close orbit 
(a = 0.0677 AU) to its host star. As shown in Table [T] its value of wqr is ~ 2°/century, while the expected value of 
Wtide is ~ 0.25 deg/century assuming fs2, P = 0.25. Unfortunately it has a high level of stellar jitter, with estimates 
ranging from 17 m sec - to 60 m sec" - 1 (jBakos et alJl2007t lLoeillet et al.ll2008| ). So even though this system has K ~ 
1000 m sec -1 it would require an unrealistically large number of observations (~ 10 4 ) per epoch in order to make its 
expected cI>gr, detectable with radial velocities. 
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The value of uj derived for HAT-P-2 b is 1.057T and we therefore would not expect to see significant variations in 
the time between primary and secondary transits in case the secondary transit was observed (see Figured]). With 
uj = 1.057T the true anomaly at the time of transit will be ft ~ 0.457T, or 1.557T. Given that the impact parameter is 
b we get from Equation [15] that \dlnD/dt\ ~ 1.67 x 10 -2 century -1 . Using the fact that D = 0.15 days we expect 
then a change of ~ 21 sec in D over 10 years due to GR. This difference could be readily detected with high precision 
photometric observations that determine D to within a, few seconds such as is possible with HST. The observations 
currently available constrain D only to within ~ 3 mins (|Bakos et al.ll2007l) , and so no first epoch suitable to measuring 
changes in D is yet in hand. 

5.2. XO-3 b 

As shown in Table [1] XO-3 b has the largest predicted wqr of all currently known exoplanets with e > 0.1. While 
its eccentricity is not as high as that of HAT-P-2 b, it is more massive and orbits closer to its host star XO-3 (also 
known as GSC 03727-01064). Asuming fc2 lP = 0.25, the expected value of wtido is ~ 11 deg/century, about three times 
as much as the contribution from GR. 

In order to assess the dctcctability of u with radial velocitie s for XO-3 b we need to know <7jitter> but unfortunately 
the precision of the radial velocity measurements presented in Uohns-Krull et all (|2008f ) is too coarse to allow a deter- 
mination o f this quantity (q b s > 100 m sec -1 ). Based on the spectral type F5V and vsini — 18.5 ± 0.2 km sec -1 
of XO-3 b (iJohns-Krull et al.ll2008f ) we can expect it to have a rather high value of stellar jitter cotter > 30 m sec -1 
(|Saar et al.lll998f) . Therefore, and just as is the case for HAT-P-2 b, we do not expect uj to be detectable with radial 
velocity observations due to the expected jitter. 

The value of uj derived for XO-3 b is consistent with and we therefore would not expect to see variations in the 
time between primary and secondary transits in case the secondary transit was observed (see Figured]). As uj ~ 0, the 
true anomaly at the time of transit will be ft ~ 7r/2 or 37r/2. Given that the impact parameter is b ss 0.8 we get from 
Equation [15] that \d\sxD/dt\ ~ 3.56 x 10 -2 century -1 . Using the fact that D = 0.14 days we expect then a change of 
~ 43 sec in D over 10 years. Just as is the case for HAT-P-2 b, this could be detectable with determinations of D to 
within a few seconds and there is no suitable first epoch yet in hand. 

6. CONCLUSIONS 

In this work we have studied the observability of the precession of periastra caused by general relativity in exoplanets. 
We additionally consider the precession caused by tidal deformations and planetary perturbers, which can produce a 
precession of comparable or greater magnitude. We consider radial velocities and transit light curve observations and 
conclude that for some methods precessions of the magnitude expected from GR will be detectable in timescales of 
~ 10 years or less for some close-in, eccentric systems. In more detail, we find that: 

1. For transiting systems, precession of periastra of the magnitude expected from GR will manifest itself through 
detectable changes in the duration of primary transit ( §3.2| or through the change in the time between primary 
and secondary transits f £)3.4p in timescales of < 10 years. The two methods are most effective at different values 
of the true anomaly. A determination of the primary transit duration time to ~ a few seconds and that of the 
secondary to ~ a minute will lead to measurable effects. The effects of GR and tidal deformations might need 
to be included in the analysis of Kepler data for eccentric, close-in systems. The transit duration of near-grazing 
systems will be particulary sensitive to changes in uj. 

2. Radial velocity observations alone would be able to detect changes in the longitude of periastron of the magnitude 
expected from GR effects only for eccentric super-massive (K ~ 1000 m sec -1 ) exoplanets orbiting close to a 
host star with a low- level of stellar jitter f H3.1|) . For the detection to be statistically significant, on the order of 
100 precise radial velocity observations are needed at each of two epochs separated by ~ 20 years. 

3. Measurements of the change over time of the period between primary transits is not currently a method that will 
lead to a detection of changes in uj of the magnitude expected from GR ( ^3.3|) . Previous works have shown that 
measuring the small difference between the radial velocity period and that of transits are not sensitive enough 
to lead to detectable changes due to GR. 

In order to contrast any detected change in the transit duration f H3.2|) or the time between primary and secondary 
(i j3.4p to the predictions of a given mechanism one needs to know the eccentricity and longitude of periastron of the 
systems, for which radial velocities are needed (although not necessarily of the precision required to directly detect 
changes in uj with them 14 ). Conversely, photometric monitoring of primary transits are useful in order to elucidate the 
nature of a detected change in uj by probing for the presence of transit time variations. The presence of the latter would 
imply that at least part of any observed changes in w could have been produced by additional planetary companions 

(ED- 

Precession of periastra caused by planetary perturbers and the effects of tidal deformations can be of comparable 
magnitude to that caused by GR f£ )2.2p . The effects of tidal deformations on the precession of periastra in particular 
may be of the same magnitude or dominate the total uj in the regime where the GR effects arc detectable ( §4.3p . While 



The eccentricity can be constrained using transit information alone, see IFord et al.l pOOl ) 
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this limits the ability to directly extract the precession due to GR given the uncertainty in the expected precession from 
tides, it might allow to study the tidally induced precession by considering the residual precession after subtracting 
the effects of GR. The latter possibility is particularly attractive in systems where the tidally induced precession 
may dominate the signal. We note that even without considering the confusing effects of tidal contributions to the 
precession, a measureme nt of ci>oR as described in this work would not be competitive in terms of precision with binary 
pulsar studies (see, e.g.. lWil ll l2006t for a review ) and would therefore not offer new tests of GR. 

The upcoming Ke pler mission expects to find a large number of massive planets transiting close to their host stars 
(jBorucki et al.ll2003f) , some of which will certainly have significant eccentricities. Furthermore, systems observed by 
Kepler will be extensively monitored for variations in their transiting time periods in order to search for terrestrial- 
mass planets using transit-time variations. We have shown that modeling of the transit time durations and further 
characterization of close-in, eccentric systems might need to take into account the effects of GR and tidal deformations 
as they will become detectable on timescales comparable to the 4-year lifetime of the mission, and certainly on follow- 
up studies after the mission ends. We have also shown that planetary companions with super-Earth masses may be 
detectable by Kepler by the change in transit durations they induce (§ 14. 2[) . Additionally, well sampled radial velocity 
curves spanning > 5 years may also be able to detect companions with super-earth massses by measuring a change in 
oj over the expected GR value for the case of super-massive, close-in systems orbiting inactive stars (§ 14.2ft . 

We would like to thank the anonymous referee for helpful suggestions and Dan Fabrycky and Andras Pal for useful 
discussions. G.B. acknowledges support provided by the National Science Foundation through grant AST-0702843. 
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